Zero-temperature equation of state of two-dimensional ^He 
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The equation of state of two-dimensional ^He at zero temperature lias been calculated using the 
diffusion Monte Carlo method. By means of a combination of the fixed-node and released-node 
techniques it is shown that backflow correlations provide a very accurate equation of state. The 
results prove unambiguously the non-self-bound character of two-dimensional '^He due to its Fermi 
statistics. We present solid evidence that the gas phase, predicted for the two-dimensional system, 
can be extrapolated to the case of ^He adsorbed on a strong substrate like graphite. 

PACS numbers: 67.55.-s, 67.70.+n, 02.70.Ss 
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■^He adsorbed on strongly interacting substrates like 
graphite, or on top of bulk ^He or *He films, constitute 
experimental realizations of quasi-two dimensional Fermi 
systems. In the last decades there has been a continued 
experimental effort to unveil the fascinating properties 
of such a nearly perfect two-dimensional Fermi liquid. 
Among these unique features, of particular relevance is 
the possibility of continuously increasing the areal den- 
sity from an almost ideal gas behavior up to a strongly 
correlated regime. This is the experimental situation 
observed for example in the two first layers of '^He ad- 
sorbed on graphite. These experimental findings indi- 
cate the non-existence of a self-bound ^He system. In 
contrast, Csathy et al. have recently studied submono- 
layer ■^He-^He mixture films on II2 and claim that '^He 
atoms appear to have condensed into a 2D self-bound 
liquid. Also, a recent theoretical study of mixture films 
points out the possibility of a dimerized ^He phase in 
a strictly 2D geometry ||]. In fact, the question of a 
self-bound 2D •^He phase has been discussed for a long 
time both from experimental (|, |5[ ^ and theoretical 
perspectives 0,110, |ll|. 

Theoretical calculations concerning the 2D ■^He system 
and the '^He films are scarce in comparison with the cor- 
responding ones for the boson isotope ''He. In addition to 
dealing with a strongly correlated system like helium, the 
Fermi statistics of "^He must be taken into account. In one 
of the pioneering works of the field, Novaco and Camp- 
bell calculated the equation of state of "^He adsorbed 
on graphite. Using lowest-order Fermi corrections, they 
concluded that the '^He film is in a gas state, contrary to 
^He which exhibits a well-established self-bound charac- 
ter [|l2, 13 1 . A comparative study of bosons and fermions 
in 2D was performed by Miller and Nosanow [|| using the 
variational method. According to their approach, and 
using a Wu-Feenberg expansion at lowest order, ^He 
cannot condense in 2D. More recently, Brami et al. ||] 
calculated the properties of a 2D ^He film using vari- 
ational Monte Carlo (VMC). They concluded that the 
presence of a transverse degree of freedom, not present 
in two dimensions, allows the system to gain enough ad- 



ditional binding energy to guarantee a liquid phase with 
a very small energy. However, a recent Creen's function 
Monte Carlo (GFMC) calculation by Whitlock et al. [|5), 
of a ''He film adsorbed on graphite, has shown that the 
energy gain with respect to the ideal 2D system is much 
smaller than the one estimated in Ref. Q. 

In this work, we use the diffusion Monte Carlo (DMC) 
method, which for bulk '^He has produced, for the first 
time, close quantitative agreement between theoretical 
and experimental results for the equation of state [ p^ . 
Our aim is to achieve the same accuracy in the present 
study of a strictly 2D ^He fluid. The energies obtained 
constitute upper bounds to the eigenvalues of the many- 
body Schrodinger equation, but the method can measure 
the quality of the bounds and provide a means of im- 
proving them. The results presented in this work show 
that backflow correlations are sufhcient to bring the sys- 
tematic error to the level of the statistical noise. The 
resulting equation of state implies the non-existence of 
self-binding, in agreement with most experimental ob- 
servations. 

DMC is a stochastic method that solves the N- 
body imaginary-time Schrodinger equation for the wave 
function f{R, t) — ip{R) 'I'(i?, t), with ip{R) a trial wave 
function used for importance sampling (see Ref. for 
a more detailed description of the actual DMC algorithm 
used) . The first and simplest approximation for the trial 
wave function tp is the Jastrow-Slater form 



(1) 



with a Jastrow factor -^j = Y[i<j exp(M(rij)) accounting 
for the dynamical correlations induced by the interatomic 
potential, and 13^ (D^) a plane- wave Slater determinant 
for the spin-up (spin-down) atoms. 

The nodal surface provided by ■0^ corresponds to the 
2D free Fermi gas. This is a first approximation since 
the real nodal surface is modified by dynamical correla- 
tions. This influence is contained in the time-dependent 
Schrodinger equation. Starting with i/j^ as zcroth or- 
der, a straightforward calculation shows that the first 
order correction to the wave function incorporates the so 
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TABLE L DMC total and partial energies (in K) of 2D ^He. 
The last column shows the upper bound (E/N)^ , relative to 
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FIG. 1: RN energies as a function of the released time t^. 
Empty and filled circles stand for ip^ and ip^^ as trial wave 
functions, respectively. Solid lines are linear fits to the calcu- 
lated data. 



called backflow correlations, a name used in analogy to 
the same type of corrections introduced by Feynman and 
Cohen in their famous paper on the ^He phonon- 
roton spectrum. In this new wave function the ar- 
guments of the orbitals (pa{i) entering and are 
shifted under the influence of the medium, 

N 

exp {ika ■ fi} = exp <^ ika ■ | n + Ab ^ v{nj)ri^ 

(2) 

The two-body correlation factor is of McMillan type, 
u{r) = — 0.5(6/r)^, and the backflow function is approxi- 
mated by a gaussian, r]{r) = exp[— ((r — rB)/i^B)^]- The 
parameters of ip have been optimized using VMC; the 
optimal values are b = 1.16cr, Ab = 0.40, re = 0.75cr, 
and LOB = 0.54(7 (cr = 2.556A). The density dependence 
of this set of parameters in the region studied here is very 
weak and can be neglected. The interatomic interaction 
corresponds to the HFD-B(HE) potential from Aziz et 
al. [|0|; its use in bulk ^^He |l§ and ^He [H has allowed 
for a very accurate calculation of their respective equa- 
tions of state. 

The upper-bound to the energy provided by the fixed- 
node (FN) approximation depends on the accuracy of the 
nodal surface oftp{R)', if it is the exact one then FN gen- 
erates the eigenvalue. Otherwise, FN yields a variational 
estimate of the energy of the system, but provides no in- 
dication of the accuracy of this estimate. In a preceding 
work on bulk '^He [ p^ we devised a combined method 
that incorporates FN, the released-node (RN) method, 
and an analytical prescription to improve ip{R). The 
RN method does not generally yield the exact energy, 
due to the growth of the boson component, but the ini- 
tial slope of the energy vs. the released time is readily 
accessible. From this slope one can guess the difference 
between the FN energy and the eigenvalue; an exact wave 
function generates zero slope. In Fig. 1, RN results at 



medium density (p — O.lOcr^^) are shown for the two 
trial wave functions tp^ and 4'^^ reported above. The 
decrease of the energy going from tp^ to ip^^ is ~ 0.01 
K, a small figure in absolute value but of the same rel- 
ative magnitude as in bulk. More importantly, the RN 
energies shown in Fig. 1 illustrate the different behavior 
of E/N with the released time t^-: the ^p^ results show a 
linear decreasing trend which disappears when backfiow 
correlations are introduced. In fact, a fit to the slope 
of E/N{ip^^) gives a value compatible with zero. This 
remarkable result can be combined with an estimation 
of the released time from E/N{ip^) which indicates that 
the systematic error of the BF energies due to the nodal 
surface is in the mK order of magnitude, i.e., it has been 
brought to the level of the typical statistical errors of this 
work. Additional insight on the high quality of the nodal 
surface provided by the inclusion of backflow correlations 
is reached by incorporating in ip the next-order analytical 
terms [Q. The results obtained with the new -0, which 
includes explicit three-body correlations in the backflow 
operator, do not show any improvements with respect to 
the BF energies. Therefore, the DMC results obtained 
for 2D •^He, using the optimized backflow correlations, 
are essentially exact; it is worth noticing that the same 
method gives unprecedent accuracy in the calculation of 
the bulk '^He equation of state. 

Results for the total and partial energies of the 2D ^He 
system as a function of the surface density are reported 
in Table I. The potential energies per particle have been 
obtained using a pure estimation method |21 in order 
to avoid any biases coming from the trial wave function. 
The kinetic energy comes from the difference between the 
total and potential energies. All the calculations have 
been carried out with 90 atoms, the finite-size simulation 
effects having been corrected for, and in practice elimi- 
nated, by summing up energy-tail contributions coming 
from both the dynamical and statistical parts. The final 
DMC energies are positive for any density and result from 
a significant cancellation between T/N and V/N , a usual 
feature in condensed helium. The last column in Table I 
contains the decrease of the energy when the nodal sur- 
face is improved by the inclusion of backfiow correlations. 
As could be expected, AE/N increases with p as does the 
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FIG. 2: Energy per particle of 2D ^He as a function of the 
density. The statistical error bars are smaller than the size 
of the symbols. The solid line corresponds to a third-degree 
polynomial fit (H) to the DMC data. 
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FIG. 3: Influence of the Fermi statistics on the energy of 
2D '^He. Filled and empty circles correspond to Fermi and 
Bose ^He, respectively. Squares represent the sum of the Bo- 
son energy and the Fermi gas kinetic energy. The lines are 
polynomial fits to data. 



relevance of correlations. 

The equation of state of 2D ^He is shown in Fig. 2 for 
the range of densities studied. In this region, our results 
are well parameterized by a cubic polynomial (solid line) 

(ii;/7V)(p) =ap + 6p2 + cpX (3) 

with optimal parameters a — 2.376(74)K(t^, h — 
-16.87(81)Kcr4, and c = 0.608(21) • lO^Kcr^. At very 
small densities (p < 0.05(7"^), the energy grows linearly 
as in a free Fermi gas {Ey/N = h'^/{2m)'Kp) but with 
different slope. At medium densities, there is a clear 
change in the slope with a flatter region suggesting the 
emergence of a minimum. Nevertheless, this minimum 
does not appear and the energy remains always positive. 
For p > 0.25cr~^ the energy increases much faster as the 
density approaches the freezing point, which experimen- 
tally is observed at p ~ 0.394o--2 g. 

The Fermi statistics of ^He atoms is the key point un- 
derlying the non-self-bound character of 2D '^He. This 
conclusion is drawn from the comparison between the 
real system and a fictitious 2D ^He boson system. In 
Fig. 3, the DMC equation of state of both systems is 
compared. Boson '^He would show a liquid phase with 
a binding energy {E/N)o = -0.1189(23) K at an equi- 
librium density po = 0.1311(17)(t~^. The lighter mass 
of the ^He atom is responsible for the large reduction of 
binding energy of boson '^He with respect to liquid *IIe: 
using the same potential, the ^He equilibrium point is 
(0. 284(7^2^ -0.897K) Some previous theoretical cal- 
culations on 2D '^He used the boson model as a reference 
system |Q, ||]. On top of this, dominant Fermi correc- 
tions were added to introduce the correct statistics. This 
perturbative approach, known as Wu-Feenberg expan- 
sion , starts at zeroth order by simply adding Ep /N to 



the boson energy. This crude estimation is plotted in Fig. 
3. As one can see, the kinetic term E-p/N leads to a large 
overestimation of E/N, and only for p < 0.03cr~^ can it 
be considered a reasonable approximation. The succes- 
sive terms in the Wu-Feenberg series show a non mono- 
tonic behavior and in general a very slow convergence: 
at medium densities {p ~ O.lcr^^) the next-order term 
is negative and approximately 100 times smaller than 
the zeroth order |jll|. Therefore, although the present 
results modify quantitatively those approximate calcula- 
tions, previous conclusions about the gas-like character 
of '^He are not altered. 

Relevant quantities from the experimental standpoint 
are the density dependence of the pressure P{p) and the 
speed of sound c{p). Both functions are shown in Fig. 4. 
They have been obtained from the polynomial fit to the 
DMC energies (^) through the thermodynamic relations 
Pip) = p^{d{E/N)/dp) and c(p) = [m-^{dP/dp)Yl^ . 
The pressure remains very low up to p ~ 0.20(t~^ and 
from then on increases much faster due to the small '^He 
mass and the rapid decrease with density of the mean 
distance between particles, due to the reduced dimen- 
sionality. An approximate estimate for the latter comes 
from 2(47rp)~i/^; at p = 0.25cr~^ it amounts only 2.88A, 
a smaller value than in bulk '^He at freezing, 4.32A. The 
speed of sound presents three different regimes as a func- 
tion of p. At very small densities, c(p) increases approxi- 
mately like p^/^ as it would correspond to a free 2D Fermi 
gas. Then, c(p) reaches a plateau up to p ~ 0.08cr^^. In 
this region, the speed of sound increases although the 
slope is very small; this behavior is a direct consequence 
of the flattening exhibited by the energy for the same 
range of densities (see Fig. 2). Finally, in the third 
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FIG. 4: Density dependence of the pressure and speed of 
sound of 2D ^He. 



be estimated from the GFMC results for the ''He hq- 
uid taking into account the different mass of the 
two isotopes in the approximate kinetic-energy correc- 
tion. For example, at densities p = 0.065cr~^, 0.131cr~^, 
and 0.170CT-2 the energy shifts are AE/N = -0.005K, 
-0.022K, and -0.037K, to be compared with the 2D en- 
ergies E/N = O.IOOK, 0.159K, and 0.220K, respectively. 
Therefore, the energy shift is by far too small to change 
the conclusion that, like the strictly 2D fluid, a thin ■^He 
film is not self-bound. 
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regime c{p) again increases with p in a more common 
way. 

The 2D system constitutes a model for a film adsorbed 
on a strong substrate like graphite. A question that nat- 
urally emerges is to what extent the presence of a trans- 
verse degree of freedom would modify the '^He proper- 
ties. First, the nodal surface could be different from 
the one of the 2D system. However, a key result of this 
work in 2D, and the previous one in 3D |16 , is that the 



real nodal surface in both cases is essentially given by 
backflow effects. In a film, the particle-particle backflow 
correlations would be mainly contained in the surface 
plane. Furthermore, both the not-in-plane contributions 
and the particle-substrate correlations, being perpendic- 
ular to the surface plane, should have a small effect on 
the backflow wave function. Therefore, the Fermi statis- 
tics of a thin '^He film can be safely considered within 
the idealized 2D geometry. Second, the additional de- 
gree of freedom perpendicular to the substrate could by 
itself lower the energy in an amount large enough to al- 
low for the existence of a liquid phase. In fact, a VMC 
calculation of '^He and ''He films adsorbed on graphite by 
Brami et al. concludes for the former the existence of 
a self-bound system with a binding energy of ^ 200mK 
at an equilibrium density ~ O.lSlu"^. To our knowledge, 
there are not DMC or GFMC calculations of ■^Hc films 
on graphite that can confirm that variational prediction. 
However, Whitlock et al. performed GFMC calcula- 
tions for ^He films adsorbed on the same substrate, and 
found a decrease in energy with respect to 2D that is 
much smaller than the results from Ref. j^. Even at the 
variational level, and using the same wave function and 
graphite-helium potential, Whitlock et al. were not 
able to reproduce the '^He energy reported in Ref. ||] 
(-0.7 vs. -1.9K ). 

The order of magnitude of the shift in energy that ap- 
pears in a ^He film with respect to the 2D system can 
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